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1. Introduction 

The past year has seen impressive efforts with respect to dynamical lattice QCD calculations 
of many hadron structure observables like form factors and moments of (generalized) parton dis- 
tribution functions, which are by now consistently carried out at pion masses as low as 300 MeV. 
New methods and techniques have been developed and successfully tested in e.g. calculations of 
the form factor and gluonic structure of the pion, exploratory studies of disconnected and strange 
quark contributions to nucleon structure observables, as well as a first investigation of the neutron 
electric polarizability with dynamical quarks. Other new and interesting developments include a 
study of recent lattice results on nucleon to A axial-vector transition form factors in the small scale 
expansion of chiral perturbation theory, an exploratory study of moments of vector meson distri- 
bution amplitudes, which represent an essential part in the understanding of rare B decays, and 
a calculation of moments of tensor generalized parton distributions of the pion, giving rise to a 
surprisingly non-trivial pion spin structure. 

Many of these observables can be defined using bi-local quark operators on the light-cone, 

Or(x) = J ^q(-^n)rW q (^n), (1.1) 

where the variable x represents the longitudinal quark momentum fraction, n is a light cone vector, 
Y = yMy 5? (jMv ? anc [ the Wilson-line % ensures gauge invariance. Lattice calculations give 
access to hadron matrix elements of x-moments of Or(x). The moments are obtained by taking the 
integral j_ x dxx n ~ l , leading to a tower of local operators 



-Mn-l 



yqfiD^---iD^- l q, (1.2) 



for the vector case r = y M , where 5? denotes symmetrization in the indices . . . and sub- 

traction of traces. Similar expressions can be obtained for the axial- vector, r = 7^/5, and tensor, 
r = <7 AIV , operators. In QCD, operators as given in Eq. (1.2) have to be renormalized and therefore 
lead in general to scale and scheme dependent quantities. Typical examples of corresponding lattice 
operators for n = 2 are discussed below in section 4.1. Non-forward matrix elements (P f \ . . . \P) of 
Eq. (1.1) can be parametrized in terms of generalized parton distributions (GPDs), e.g. H{x,^,t) 
and £(jc, t, , t) for the nucleon in the vector case, 

(P f \0^x)\P) = U(P f )^fH(x,^t) + ^^E(x,^t)jU(P), (1.3) 

where t — A 2 = (P 1 — P) 2 is the momentum transfer squared, and £ corresponds to the longitudinal 
momentum transferred to the hadron (for reviews on GPDs, see [1, 2]). Note that H{x 1 1; , t) reduces 
in the forward limit, A — > 0, to the usual unpolarized parton distribution function, H(x, 0, 0) = q(x). 
Corresponding nucleon matrix elements of the local operators in Eq. (1.2) for n = 1,2 are given by 



(P'\^\P) = U(P')I^F l ( t ) + ^—tF 2 (t))U(P), 



(p'\^\p) = yu(P')[r^A 20 (t) + ^^P^B 20 (t) + ^-C 20 (t) 
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Figure 1: Results for gA from QCDSF/UKQCD [6]. Symbols are explained in the text. 

where P = (P' + P)/2, and are the Dirac and Pauli form factors. Since ffy^ 1 corresponds to 
the quark contribution to the QCD energy momentum tensor, the generalized form factors (GFFs) 
A2o(0 an( i ^20 (0 at f = can be directly related to the quark momentum fraction, (x) q = A2o(0), 
and the total quark angular momentum, J q = (A2o(0) + 2?2o(0))/2. The lattice calculation of matrix 
elements as in Eq. (1.4) for the nucleon from two- and three-point functions, and the subsequent 
extraction of moments of parton distributions and (generalized) form factors follows standard meth- 
ods described in e.g. [3, 4, 5] and references therein. Some new techniques and strategies used in 
the calculation of hadron structure observables on the lattice will be discussed in the sections below. 

2. Form factors 

2.1 Axial- vector coupling constant 

The axial- vector coupling constant ("axial charge") of the nucleon, gA = 1.2695(29), is ex- 
perimentally very well known from neutron beta decay. It corresponds to the forward limit of 
the isovector axial-vector form factor, Ga(Q 2 — > 0) = gA- In a lattice calculation, using isospin 
symmetry, gA can be extracted from the proton forward matrix element of the u — d axial-vector 
current, 

(p,s\uY^r 5 u-dr^r 5 d\p,s) = u(p,s)^r 5 u(p,s) gA . (2.1) 

Equation (2. 1) shows that g A is directly related to the isovector quark spin fraction, g A = AZ M _j, and 
that disconnected diagrams, which cancel out for the flavor combination u — d, do not contribute. 
The chiral limit value of the axial- vector coupling, g° A , is a fundamental low energy constant of 
the chiral effective field theory of QCD (ChEFT). Based on the Adler-Weisberger sum rule, one 
might anticipate that the A resonance, in addition to pion and nucleon degrees of freedom, plays an 
important role in the low energy description of gA in the framework of chiral perturbation theory 
(ChPT). It is also known that the axial- vector coupling is particularly sensitive to finite volume 
effects in, e.g., a lattice simulation. The application of ChPT including the A resonance in a fi- 
nite volume to the description and extrapolation of lattice data on g A therefore seems to be very 
promising. Figure 1 shows new preliminary results from QCDSF/UKQCD for gA, based on Nf = 2 
improved Wilson fermions and Wilson glue, for pion masses as low as ~ 340 MeV [6]. The opera- 
tor has been non-perturbatively renormalized using the Rome-Southampton method. The scattered 
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open symbols on the left already indicate that finite size effects (FSEs) may be large. Indeed, a 
fit based on ChPT including explicitly the A resonance in a finite volume [7, 8, 9, 10] reveals a 
significant dependence on the lattice extent, L, see the RHS of Fig. 1 for a fixed pion mass of « 600 
MeV, and allows for a simultaneous description of the m n and L dependence of the lattice data. 
The lattice results and the corresponding chiral fit, 
projected to the infinite volume, are represented 
by the filled symbols and the shaded error band, 
respectively, on the LHS of Fig. 1. The curvature 
of the chiral extrapolation is just strong enough 
as to provide an overlap with the experimental 
value at the physical point, indicated by the star. 
Preliminary results on g A by RBC-UKQCD [11] 
based on Nf = 2+1 domain wall fermions with 
L s = 16 and volumes of 1.8 and 2.7 fm are shown 
in Fig. 2. It is gratifying to see that dynamical chi- 
ral fermion calculations at pion masses as low as 
c± 330 MeV lead to results for gA with statistical 
errors around the 5% level. However, the signif- 
icantly lower result at the lowest pion mass, rep- 
resented by the shaded circle, may indicate that FSEs are substantial. For comparison with results 
from LHPC at similar pion masses see [12]. 



m'[GeV 2 ] 



Figure 2: Results for gA from RBC- 
UKQCD[11]. 



2.2 Nucleon charge radii and anomalous magnetic moments 

Despite interesting developments on the theoretical as well as the experimental side in this 
field, publications on nucleon form factors (FFs) F\^{Q 2 ) in unquenched lattice QCD are surpris- 
ingly scarce (for recent reviews, see [13, 14]). This includes such fundamental observables as mean 
square charge radii, (r 2 ), and the anomalous magnetic moment, K*, 



(r 2 )i 



d 



Fi(Q 



J<2 2 =o 

K = F 2 (Q 2 = 0) = ! 1-F 1 (Q 2 = 0), 



(2.2) 
(2.3) 



where / = 1,2 and ji = Gm{Q 2 = 0) is the nucleon magnetic moment. Preliminary results on 
(r 2 ) u x ~ d as a function of the pion mass by QCDSF/UKQCD [15] are shown in Fig. 3. The mean 
square isovector charge radius has been obtained from a new parametrization [15] of the Q 2 - 
dependence of the lattice results for the Dirac FF F\ (Q 2 ). Even at the lowest pion mass of ~ 340 
MeV, the lattice data points are almost a factor of two below the experimental result, which is rep- 
resented by the star. Chiral perturbation theory in the form of the small scale expansion (SSE) to 
<^(e 3 ) in the infinite volume [16] predicts a rather strong slope in rn n , which barely connects the 
leftmost lattice datapoints with the experimental number, as indicated by the dashed line. It will be 
highly interesting to see if the lattice results follow the chiral extrapolation curve at pion masses 
below 300 MeV. 

Similar results have been obtained by RBC-UKQCD [11], Fig. 4, where the charge radius has 
been extracted from a dipole fit to the lattice results, i.e. (r 2 ) u l ~ d = \2jm 2 D , with the dipole mass 



4 



Progress in hadron structure physics on the lattice 



Philipp Hagler 



0.5 



500 1000 
ni / MeV 



Figure 3: Isovector mean square radius (r 2 )i from 
QCDSF/UKQCD [15]. 
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Figure 4: Isovector root mean square radius {r ) x 
from RBC-UKQCD [11]. 





Figure 5: Quark contributions to the nucleon anomalous magnetic moment k from QCDSF/UKQCD [15]. 



mo- Figure 5 shows preliminary results from QCDSF/UKQCD [15] on the up- and down-quark 
contributions to the nucleon anomalous magnetic moment. The results of 3-parameter chiral fits 
based on the small scale expansion (SSE) [16] to the lattice datapoints are represented by the shaded 
errors bands, showing good agreement with the experimental values, indicated by the stars, at the 
physical point. Nevertheless, these results should be taken with due caution since contributions 
from disconnected diagrams have not been included. 

Further results on vector and axial-vector nucleon form factors, including the axial-vector 
charge radius and the pseudoscalar coupling constant, have been presented by QCDSF/UKQCD 
[15], RBC-UKQCD [11] and the Athens-Cyprus-MIT collaboration [17, 18]. 

2.3 Pion form factor 

Substantial progress based on new methods and techniques can be seen in lattice calculations 
of, e.g., the pion form factor F K (Q 2 ). The analysis of many important hadron structure observables 
requires small but non-zero values for the (squared) momentum transfer q 2 = (pf — pi) 2 = —Q 2 . 
Typical examples are charge radii (2.2), anomalous magnetic moments (2.3) and angular momen- 
tum contributions to the nucleon spin (see section 5.1), where in the latter two cases the corre- 
sponding (generalized) FFs cannot be obtained at Q 2 = because their contribution to hadron 
matrix elements vanishes in the limit Q 2 — ► 0. This poses a problem for current lattice calcu- 
lations, since due to limited spatial lattice extents, L = 16, ... ,48, the lowest available non-zero 
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Figure 6: Pion FF using pTBCs from [21]. 
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Figure 7: Pion FF from ETMC [22]. 



lattice momenta are p = 2n/(aL) « 300 MeV. A way around this in a dynamical lattice calcu- 
lation is to use partially twisted boundary conditions (pTBCs) [19, 20], leading to a modified 
expression for the momentum transfer, q\ BC = (E^tbc ~ Ei,tbc) 2 — (? + (0/ — 6i)/{ a L)) 2 with 
Ei(f),TBC = ((Pi(f) + ®i(f)/( a L)) 2 + m 2 y/ 2 , where the twisting angles (0/(/))j=i,2,3 can be tuned 
continuously. An exploratory study of this technique, based on Nf = 2 + 1 DW fermions, has been 
published recently [21], and the results for the pion FF are shown in Fig. 6. The vertical dashed 
lines indicate the lowest non-zero q 2 which can be obtained in the conventional calculation, and all 
lattice datapoints to the left of these lines (except at q 2 = 0) could only be obtained using pTBCs. 
The smallness of the error bars for non-zero twisting angle but vanishing initial and final hadron 
three-momenta (two leftmost filled squares) is impressive, in particular in comparison with the er- 
ror bars of datapoints obtained for non-zero hadron momenta, which are a factor of 5 to 10 larger 
due to larger statistical noise. Similarly noteworthy progress has been reported by ETMC in the 
calculation of F n (Q 2 ) based on twisted mass Wilson fermions [22] at pion masses as low as 300 
MeV. All-to-all propagators have been evaluated using stochastic sources and utilized for the calcu- 
lation of pion three- and two-point-functions, leading to a substantial improvement of the precision 
compared to the standard approach employing point-to-all propagators, see Fig. 7. 

A calculation of the pion FF at pion masses 
down to 288 MeV, using Nf = 2 overlap fermions 
with a Wilson kernel and the Iwasaki gauge ac- 
tion in a fixed topological sector, has been pre- 
sented by JLQCD [23]. The calculation of the 
pion correlation functions is based on all-to-all 
propagators, which have been evaluated following 
the strategy proposed in [24]. Preliminary results 
for the pion charge radius, obtained from a vari- 
ety of parametrizations of the Q 2 -dependence of 
Fti(Q 2 ), versus the pion mass are shown in Fig. 8. 
It is remarkable that a statistical precision at the 
few percent level has been achieved with chiral 
fermions in a dynamical calculation at such low pion masses. Further studies of systematic effects 
related to e.g. the finite lattice volume, setting of the scale and fixed topology may help to under- 
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Figure 9: The transition FFs Cj 6 calculated in [17]. Figure 10: Cj in ChPT (SSE) by Procura [28]. 



stand why the absolute values for {r 2 ) K in Fig. 8 are rather low compared to experiment, and also 
lower than the results obtained recently in an extensive study of F n (Q 2 ) based on Nf = 2 improved 
Wilson fermions [25]. 

2.4 Vector meson quadrupole moments 

Matrix elements of the vector current for spin one hadrons, e.g. (p(P' 1 S f )\qy !1 q\p(P^S)} for 
the p, can be parametrized by the three form factors Ge(Q 2 ), Gm(Q 2 ) and Gq{Q 2 ). Their forward 
values at Q 2 = are directly related to the electric charge q, the magnetic moment /i, and the 
quadrupole moment Q, respectively, e.g. Gq(Q 2 = 0) = m 2 Q p . All these observables obviously 
yield important information about hadron structure, and quadrupole moments in particular can 
exhibit possible spatial deformations of hadrons. Recently, the vector meson form factors have been 
studied in quenched lattice QCD using fat-link irrelevant clover (FLIC) fermions at pion masses 
as low as « 290 MeV [26]. As in the case of magnetic moments, quadrupole moments cannot 
be directly extracted from the lattice correlators, since their contribution to the matrix elements 
vanishes in the limit Q 2 — ► 0. Assuming that the <2 2 -dependence f the electric and quadrupole 
form factors is similar at low Q 2 , Gq(Q 2 )/Ge(Q 2 ) ~ const., one obtains an approximation for the 
quadrupole moment of, e.g., the p from Q p « Gq{Q 2 ) / (G^(<2 2 )m 2 ) for small Q 2 . Based on results 
for G E (Q 2 ) and G Q (Q 2 ) at Q 2 ~ 0.22 GeV 2 , a clearly non-zero small negative value, Q p « -0.007 
fm 2 , has been obtained for the p quadrupole moment [26]. Together with results for the charge 
radius (r 2 ) p of the p, this leads to a ratio of \Q p \/(r 2 ) p w 1/30 — 1/50, indicating that the p is not 
completely spherically symmetric but slightly oblate. 

2.5 Nucleon to A axial-vector transition form factors 

In contrast to the vector mesons discussed in the previous section, the nucleon as a spin 1/2 
particle has no (static) quadrupole moment which could be utilized to measure possible deviations 
from spherical symmetry. Instead, spin 1/2 to spin 3/2, in particular nucleon to A, vector transi- 
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tions can and have been used to study possible non-zero quadrupole amplitudes, see, e.g., [27] and 
references therein. Recently, the corresponding axial-vector nucleon to A transition form factors 
have been investigated in detail in lattice simulations based on Wilson fermions in the quenched and 
unquenched theory and in a hybrid approach of Nf = 2 + 1 DW valence quarks on top of AsqTad 
staggered sea quarks (MILC configurations) [17, 18]. The results for the two dominant transition 
form factors, C^(Q 2 ) and C£(£ 2 ), are shown in Fig. 9. Apart from the overall very good statistical 
precision and agreement between the different lattice actions, three lattice data points at low Q 2 and 
the lowest pion mass ~ 360 MeV (crossed squares) stand out in the sense that they do not follow 
the general trend of the data. In addition to studying possible systematic uncertainties of the cal- 
culation, it is certainly very interesting to see if this downwards bending at low Q 2 and m 2 can be 
understood within ChPT. This has been looked at recently by Procura [28] in the framework of the 
small scale expansion (SSE) to leading one-loop accuracy ff{p 3 ). The result for C$ is, in short hand 
form, given by C^' SSE (2 2 ,m^) = cf' + C^m* +^ Q Q 2 + ^ loop (m 7l ,c A ,g A ,g l ,/ w ,A), show- 
ing to this order a linear dependence on Q 2 , while the non- analytic dependence on m K is hidden in 
C^J° op . The low energy constants in ^ loop ma y be fixed to known values from the literature and 
by a fit to lattice data for the axial- vector coupling [9]. The result of a three parameter fit (with 
parameter C^'° plus the two counter terms °c Q 2 and oc m 2 ) to six lattice data points of Fig. 9 is 
presented in Fig. 10. It is encouraging to see that, within the error bands indicated by the dashed 
lines, the SSE calculation is compatible with the results from the hybrid lattice calculation. 

3. Polarizabilities in a dynamical calculation 

Nucleon electric and magnetic polarizabilities, <Xe and /3m, parametrize the forward Compton 
scattering amplitude at ff{(0 2 ) in a low energy expansion in the photon energy CO. They describe 
the resistance of the internal degrees of freedom of the nucleon against external forces created by 
the electromagnetic field of the scattered photon and therefore encode significant information on 
nucleon structure. The polarizabilities are related to a shift of the hadron mass through the effective 
Hamiltonian H eff = -(a E E 2 + j8 M £ 2 ) /2. 
In a recent study, Engelhardt [29, 30] ex- 
plored the feasibility of a new method to 
compute the electric polarizability of the 
neutron in a dynamical lattice calculation. 
The inclusion of external fields in the calcu- 
lation of gauge configurations renders such 
dynamical calculations at first sight pro- 
hibitively expensive. A way around this is 
to expand the neutron two-point function in 

the external (electric) field E to the desired order, &(E 2 ), resulting in specific four-point functions. 
Examples for connected and disconnected contributions are given in Fig. 11, where a cross repre- 
sents an insertion linear in the external field. The mass shift, and thereby the electric polarizability, 
can then be extracted from the Aslope of the ratio of the expanded two-point function to the two- 
point function for vanishing external fields, Ri{t) = G^ (t) / G(£= ) (0' where / is the source-sink 
separation in Euclidean time direction. The numerical calculation has been based on the hybrid 



o 



■ + — I o- K) |+ — ! cH ho ! + •■ 

2 ' i i ' 2 \ \ 

W 102(12) W W J01(2) W 



Figure 11: Examples for connected (left) and discon- 
nected (right) four-point function contributions to the 
neutron electric polarizability [29] . 
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Figure 12: Contribution to the ratio R2 from dia- 
gram J01 (see Fig. 11) [29]. 



Figure 13: Contribution to the ratio R2 from dia- 
gram 102 (see Fig. 11) [29]. 



approach described in [12, 31], i.e. domain wall valence fermions and a MILC gauge ensemble. 
Examples for ratios are shown in Fig. 12 and Fig. 13 for a pion mass of ~ 760 MeV. The slope 
of R2 for the disconnected contributions, which have been evaluated using stochastic sources, is 
compatible with zero in Fig. 12, while non-zero slopes can be observed in Fig. 13 for individual 
connected contributions. In contrast to previous works, great care has been taken to disentangle 
the polarizability from additional effects due to constant electromagnetic fields on a finite lattice, 
resulting in = —2.0(9) • 10~ 4 fm 3 for the electric polarizability of the neutron at m % c± 760 MeV. 
The substantial difference to the experimental value = 11.6(1.5) • 10~ 4 fm 3 may be explained 
in the framework of chiral perturbation theory, which predicts a strong pion mass dependence of 
the form \jm % at low m % . For comparison, we refer to a previous investigation in the quenched 
approximation in Ref. [32]. 

4. Moments of parton distribution functions 

4.1 Quark momentum fractions in the nucleon 

The momentum fractions carried by the quarks and gluons in the nucleon, (x) qig , represent 
fundamental hadron structure observables. They are directly related to the moments of structure 
functions through the operator product expansion (OPE). Although separately dependent on the 
scale, /I, and the renormalization scheme, they must add up to one according to the fundamental 
momentum sum rule, Y, q ( x )q + ( x )g = 1» independent of the scheme and scale. The quark momen- 
tum fractions in particular are very well known most notably from deep inelastic scattering (DIS) 
experiments in combination with phenomenology, and have been studied already for a long time in 
lattice QCD. 

As discussed in the introduction, leading twist continuum operators, Eq. (1.2), are used to 
define moments of (generalized) PDFs. On a discrete space-time lattice, the corresponding traceless 
and (anti-) symmetrized operators are classified according to irreducible representations of the 
hypercubic group H(4). Typical lattice operators related to the quark momentum fractions are 
&w2a = ^7{i^>4}^ an d &v2b = #(7{4^>4} l/^ELi 7{iA})#> belonging to the representations t| and 
T 3 , respectively. The corresponding lattice quark momentum fractions are often denoted by \2a and 
\2b> In the continuum limit, both operators should of course yield the same result. Figure 14 shows 
preliminary results from QCDSF/UKQCD [6] for the isovector momentum fraction v^ versus 
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Figure 14: Quark momentum frac- Figure 15: Quark momentum fraction (x) M _j in the nucleon 

tion {x) u -d in the nucleon from from LHPC [31, 35]. 

QCDSF/UKQCD[6]. 



the pion mass squared at a fixed coupling j8 = 5.29. The operator has been non-perturbatively 
renormalized following the Rome-Southampton method and transformed to the MS scheme at a 
scale of 4 GeV 2 . Since the dependence on the lattice spacing appears to be small [6], we may 
directly compare with results from phenomenology, e.g. (x)^ 1 ^ 6 ~ 0.155, indicated by the star. 
At the lowest pion mass of ~ 340 MeV, the lattice data point is « 40% above the CTEQ value. In 
order to see wether this gap can be bridged, recent results from covariant baryon ChPT (CBChPT) 
[33, 34] have been used to fit the lattice data. One advantage of the covariant formalism is that all 
powers of the ratio {m n jm^) are resummed, in contrast to the more widely used non-relativistic 
heavy baryon approach, which is based on a simultaneous expansion in m K /A x and m^/m^, where 
A x ~ 1 GeV is the chiral symmetry breaking scale. The result of a 2-parameter CBChPT-fit to the 
lattice data is shown as shaded error band in Fig. 14. Although a rather strong downwards bending 
is visible, the fit misses the phenomenological value at the physical point by « 20%. 

Results for the isovector momentum fraction by LHPC [31, 35], obtained in the framework 
of an extensive study of moments of GPDs, are presented in Fig. 15. Here it is useful to note 
that the momentum fraction is equal to the (n — 2) -moment of the vector GPD H in the forward 
limit, / dxxH{x 1 Q 1 t = 0) = A2o(t = 0) = (x). The calculation is based on a combined analysis of 
the perturbatively renormalized (with NP-improvement, see [31]) lattice operators in the represen- 
tations t| and %\ . The shaded error band is the result of a global CBChPT- fit to the lattice data 
for m 2 < 0.47 GeV 2 and momentum transfers \t\ < 0.47 GeV 2 . A very good agreement with the 
phenomenological value is found at the physical pion mass. This may be attributed to the overall 
lower normalization of the lattice data points, in particular compared to the results in Fig. 14, as 
well as the inclusion of all powers {m n /m^) n in the CBChPT approach, compared to the heavy 
baryon formalism, as can be seen from the dotted and dashed lines in Fig. 15. Even taking into ac- 
count that the calculations presented in [6] (Fig. 14) and [31] (Fig. 15) are based on different lattice 
actions and methods, the difference in the normalization of (x) u -d can be regarded as substantial 
and certainly demands a closer examination, including possible finite size effects [36]. 

New results on isovector momentum fractions in the nucleon based on dynamical DW fermions 
have been presented by RBC-UKQCD [11]. 
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Figure 17: Gluon momentum fraction (x)g in a "heavy" 
pion with m K ~ 1.06 GeV [44]. 



In addition to isovector quantities, which provide essential information on hadron structure, 
calculations in the isosinglet and the strange-quark sector are required for the analysis of a large 
number of important observables and e.g. the momentum and spin sum rules. However, the numer- 
ically expensive contributions from disconnected diagrams are mostly neglected so far, and such 
calculations are therefore subject to systematic uncertainties of largely unknown size. Preliminary 
results on the disconnected contributions to the quark momentum fraction in the nucleon have been 
obtained in a recent effort by Deka and Liu [37], using the standard Wilson gauge and fermion 
action in the quenched approximation for pion masses down to 480 MeV. Their calculation of the 
all-to-all propagators is based on the Z2 noise estimator and an unbiased subtraction to reduce the 
variance using the hopping parameter expansion. Employing perturbatively renormalized opera- 
tors of the off-diagonal type, e.g. & Y 2a, they find in the MS scheme at a scale of ~ 3 GeV 2 a 
strange quark contribution to the momentum fraction of (x) s = 0.046(17), and (x)^f = 0.055(18) 
for the contribution from disconnected insertions. These results have been linearly extrapolated 
to the chiral limit. This shows that disconnected diagrams may give a substantial, non-negligible 
contribution to the total quark momentum fraction as large as ~ 10% of the value for the connected 
part in the quenched approximation, (x) c ™ d « 0.65. Studies along similar lines have been been pre- 
sented by Babich [38], using stochastic noise estimators to compute the strange quark contribution 
to the axial- vector and scalar form factor of the nucleon, and Collins [39], discussing the use of 
stochastic methods combined with unbiased subtraction and a truncated inversion of the stochastic 
propagator in the calculation of disconnected contributions to hadronic structure. 

A possible alternative to the calculation of moments of PDFs using local operators is the direct 
study of the OPE on the lattice. Interesting preliminary results in this direction based on overlap 
fermions have been presented by QCDSF/UKQCD [40]. 

4.2 Parton momentum fractions in the pion 

Remarkable work in progress by ETMC [41] on the momentum fraction of quarks in the pion, 
i.e. (x)% + , is presented in Fig. 16. The calculation is based on the if (4) -operator ff^b defined in 
the previous section, and the pion 2- and 3-point functions have been evaluated using stochastic 
sources instead of point sources. Using just ~ 300 configurations for the measurement, this still 
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leads to a high numerical precision at the lowest pion mass of ~ 300 MeV (cross and leftmost 
square in Fig. 16), where the statistical error for the unrenormalized quark momentum fraction in 
the pion is <^(3%). These results, as well as calculations by QCDSF/UKQCD on the pion structure 
[42], indicate that the quarks in the pion carry around 50 — 60% of the total pion momentum. 

This may serve as a reminder of the well-known fact that gluons give substantial, if not major, 
contributions to such fundamental properties as the mass and momentum of hadrons. The calcu- 
lation of gluonic contributions in lattice QCD to hadron structure observables, however, poses a 
long-standing problem, since short-range quantum fluctuations lead to noisy signals. Using a large 
number of configurations in a quenched Wilson fermion and Wilson glue simulation, this problem 
was tackled in a pioneering exploratory study of the gluon momentum fraction in the nucleon in 
Ref. [43] more than ten years ago. Recently, Meyer and Negele [44, 45] have studied the gluon 
momentum fraction in the pion, (x)^, in quenched lattice QCD using the Wilson action. In order 
to reduce the gauge-field fluctuations, they employed HYP smearing to smooth the fields. By com- 
paring bare-clover, HYP-plaquette and HYP-clover discretizations of the gluon energy momentum 
tensor, they found that the HYP-plaquette version reduces the variance of the entropy density by a 
factor of ~ 40 relative to the bare-plaquette discretization, without leading to uncontrollable non- 
locality effects. Their results for the bare gluon momentum fraction are shown in Fig. 17 as a 
function of the minimum distance xo in Euclidean time direction between the gluon operator and 
the pion source and sink. Taking into account the mixing of the singlet quark and gluon operators 
in the renormalization procedure, they obtain (x)g = 0.37 ± 0.08 stat ± 0. 12 sys in the MS scheme at a 
scale of 4 GeV 2 and a pion mass of m K ~ 0.89 GeV. The estimated systematic error of ±0.12 orig- 
inates from the fermion normalization constant Zf which is so far not known beyond the (trivial) 
tree level value, Zf = l + &{g 2 ). Together with corresponding results from [46] for the connected 
quark momentum fraction in the pion, they find Y,q( x )q ,C ° n + ( x )g = 0.99 ± 0.08 stat ± 0.12 sys , very 
close to one, but still allowing for disconnected quark contributions as large as 10% within errors. 

4.3 Moments of meson distribution amplitudes 

Important information on hadron structure for small transverse parton separations is provided 
by hadronic distribution amplitudes (DAs). Meson DAs are denoted by 0(£), where t; is directly 
related to the quark and anti-quark longitudinal momentum fractions through x= 1/2(1 + ^) and 
l—x= 1/2(1 — £), respectively. Following calculations of the lowest two moments of pseudo- 
scalar meson DAs [47], QCDSF/UKQCD now presented first results for moments of vector meson 
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Figure 20: Decomposition of the nucleon spin 
from QCDSF/UKQCD [50]. 



DAs, specifically the K*, in lattice QCD [48]. The DA (j>K*(%) is an essential ingredient in the 
description of rare B decays like B —> K*l + l~ in the framework of QCD factorization. Since these 
decays are induced by flavor-changing neutral current transitions only occurring through loops in 
the standard model, they are an ideal tool to probe new physics. Moments of meson DAs, (£, n ~ l ) = 
J parametrize meson-to-vacuum matrix elements of the quark operators m 
Eq. (1.2). For the calculation of the (n = 2)-moment (%)k*, the (NP-renormalized) operators G^a 
and & V 2b, discussed in section 4.1, have been used. Figure 18 shows the result for (%) a ,K* at a fixed 
sea quark mass and coupling versus the difference m\ — m\ in the squared pseudo-scalar lattice 
masses. The vertical dashed-dotted line represents the physical value for m\ — m\. An interpolation 
of (Z)a,K* 1° this point, followed by a linear chiral extrapolation in m|, gives (%) a ,K* ^0.036(3), 
in the MS scheme at a scale of 4 GeV 2 . This first- time study proves the feasibility of future high 
statistics dynamical lattice QCD calculations of moments of vector meson DAs. 

Promising results on moments of pion and kaon DAs in a dynamical DW lattice calculation 
with pion masses as low as 300 MeV have been presented by Sachrajda [49] for RBC-UKQCD. 
Using perturbatively renormalized operators, the preliminary results in the MS scheme for fi = 2 
GeV are {^ 2 ) n « 0.28(3) for the pion, and (%) K w 0.029(2), (% 2 ) K « 0.27(2) for the kaon, showing 
that "lattice calculations can be performed with an excellent precision" [49]. 



5. Moments of generalized parton distributions 

A number of recent lattice results on moments of GPDs have been presented by QCDSF/UKQCD 
for the nucleon [50] and the pion [42] and by LHPC for the nucleon [35]. Below, we discuss only 
a small selection lattice of these results, related to the nucleon and pion spin structure. 

5.1 Decomposition of the nucleon spin 

Of particular interest in the investigation of the structure of the nucleon is the decomposition 
of the nucleon spin 1/2 in terms of quark and gluon spin and orbital angular momentum (OAM) 
contributions. According to the nucleon spin sum rule , we have 

± = lAL + AG + L q + L gl (5.1) 
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where AZ and AG are the standard gauge-invariant quark and gluon spin fractions, while the OAM 
are defined by L q = J q — AZ/2 and L q = J g — AG. The gauge-invariant total angular momenta J qig 
are related to the forward values of the (n = 2) -moments of the GPDs H and E through 2J qig = 
J dxx(H(x^,0) J r E(x^,0)) = (A|(f (0) + B%§ (0)) [51]. It is important to note that the decompo- 
sition of the nucleon spin is in general scale and scheme dependent. Since the A^f (/ = 0) = {x) qig 
are just the well-known quark and gluon momentum fractions, the genuine contribution from GPDs 
to the total nucleon spin is given by the GFFs fi|(f ■ Figure 19 shows preliminary results from 
QCDSF/UKQCD for the isovector B u 2 ~ d {t = 0) as a function of the pion mass squared [50]. The 
lattice data points have been chirally extrapolated based on results from the covariant BChPT cal- 
culation presented in [33, 34], as indicated by the shaded band. A sizeable contribution of B^ d {t = 
0) = 0.269(20) is found at the physical point. Corresponding results in the isosinglet channel, to- 
gether with preliminary results for AL u+d and {x) u+ d lead to the nucleon spin decomposition shown 
in Fig. 20, in the MS scheme at a scale of 4 GeV 2 . It is remarkable that the u + d-quark OAM 
contribution is compatible with zero over the 
full range of accessible pion masses. A fur- 
ther decomposition in terms of up and down 
quark spin and OAM contributions to the nu- 
cleon spin by LHPC [31, 35] is presented in 
Fig. 21. Interestingly, the separate u and d 
quark OAM contributions are quite sizeable 
but opposite in sign, Lj ~ — L u « 30% of 
1/2, and only cancel in sum, L u+ d ~ 0, for 
all pion masses. Similarly, the down quark 
spin and OAM contributions in Fig. 21 are 
of same size but opposite in sign, so that the 
total down quark angular momentum is ap- 
proximately zero, Jd = AZj/2 + Lj « 0. The 
chiral extrapolations in [31] indicate that this 
also holds at the physical point, so that the total quark contribution to the spin of the nucleon is to 
nearly 100% coming from the up quarks. All these results should be taken with due caution, since 
disconnected contributions have not been included in the calculations presented in Fig. 20 and 
Fig. 21. A comparison with recent combined experimental results on J u ^ from JLab and HERMES 
can be found in [53]. 

5.2 Spin structure of the pion 

A longstanding, important question is how quarks are spatially distributed inside hadrons. In 
[54], it has been shown that vector and tensor GPDs in particular allow for a proper definition of 
transverse coordinate (impact parameter) space densities of transversely polarized quarks in the nu- 
cleon. Impact parameter densities, which can be determined from Fourier transforms of GPDs with 
respect to the transverse momentum transfer, H(x, b 2 ^) = (2n)~ 2 J d 2 A± exp(-iAj_fcj_)//(jt, 0, -Aj_), 
have been introduced by Burkardt [55] and are illustrated in Fig. 22. In analogy to the nucleon case 
[56], the density p(x,b±,s±) of quarks with longitudinal momentum fraction x, transverse spin s± 
at transverse distance b± of the center of momentum in the pion can be directly obtained from the 
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Figure 21: Decomposition of the nucleon spin from 
LHPC [31]. Filled stars represent recent HERMES 
results for AL u ' d /2 [52]. 
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Figure 22: Density interpretation of Figure 23: Densities of up-quarks in the 7T + from 
GPDs in impact parameter space. QCDSF/UKQCD [42, 57]. The arrow indicates the orientation 

of the transverse quark spin s± . 

pion vector and tensor GPDs, H^^x^b^) and E^x^b]), respectively. The lowest two moments of 
these pion GPDs have been calculated recently by QCDSF/UKQCD [42, 57], and the correspond- 
ing results for the lowest x-moment of the density, p n=l (b±,s±), are shown in Fig. 23 for up quarks 
in a 7T + . Compared to the unpolarized case on the left, the density of quarks with transverse spin 
in x-direction on the right in Fig. 23 is strongly distorted in ^-direction. This result is very similar 
to what has been obtained for quarks in the nucleon [56] and proves the presence of strong correla- 
tions of the form s± x b± between transverse spin and coordinate degrees of freedom of quarks in 
the pion. Following arguments by Burkardt [58] and others [59], this distortion also indicates that 
the so-called Boer-Mulders function ti^' (x,k±) of the pion, describing correlations of the intrinsic 
quark transverse momentum, k±, and the transverse quark spin, is large and negative. This may be 
important for future studies of azimuthal asymmetries in polarized and unpolarized Tip Drell-Yan 
production at COMPASS/CERN. 

Preliminary results of a first direct study of transverse momentum dependent PDFs on the 
lattice based on non-local operators have been presented by Musch [60]. 

6. Summary 

Recent lattice hadron structure calculations have provided substantial new insights into the 
shape of hadrons, their momentum and spin structure in terms of quarks and gluons and the spatial 
distribution of quarks in the nucleon and the pion. As a consequence, an increasing number of lat- 
tice results in this field have a direct impact on current and future experimental measurements and 
related phenomenology. Remarkable progress has been made in calculations of the pion structure 
based on dynamical chiral fermions, using all-to-all propagators as well as partially twisted bound- 
ary conditions. It will be interesting to see if this can be repeated in the near future for some of the 
more involved nucleon structure observables. Despite these successes, several questions remain to 
be answered, regarding, e.g., discrepancies between different computations of the isovector quark 
momentum fraction in the nucleon, the consistently low lattice results for mean square charge radii, 
and the importance of disconnected contributions in dynamical simulations. Finally, as we enter the 
regime of precision calculations at pion masses ~ 300 MeV and below, chiral extrapolations based 
on ChEFT of QCD become increasingly relevant and will help to significantly reduce systematic 
uncertainties of the lattice computations. 
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